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Abstract 

Centre manifold theory is applied to some dynamical systems arising 
from spatially homogeneous cosmological models. Detailed information is 
obtained concerning the late-time behaviour of solutions of the Einstein 
equations of Bianchi type III with collisionless matter. In addition some 
statements in the literature on solutions of the Einstein equations coupled 
to a massive scalar held are proved rigorously. 


1 Introduction 

Over the years a great deal of effort has been put into the study of homogeneous 
solutions of the Einstein equations coupled to various matter sources. This class 
of solutions has the advantage that the evolution equations reduce to ordinary 
differential equations. This allows a mathematically rigorous treatment of many 
issues for which corresponding results seem out of reach at present when partial 
differential equations are involved. The area of mathematics which is relevant 
for the study of the time evolution of homogeneous cosmological models is the 
theory of dynamical systems. An excellent source of background information 
on the application of dynamical systems to cosmology is the book edited by 
Wainwright and Ellis |]l3| . 

When a dynamical system (i.e. a system of ordinary differential equations) 
constituting a model of an aspect of the real world has been derived a natural 
task is to understand the qualitative behaviour of general solutions of the equa¬ 
tions. Doing this may involve combining many techniques. In the following we 
concentrate on one particular type of technique. In understanding a dynami¬ 
cal system it is useful to determine the stationary solutions. Once these have 
been found it is desirable to know the behaviour of solutions which remain in a 
neighbourhood of one of these stationary solutions. If the system is linearized 
about the stationary point then there is often a relationship between solutions 
of the linearized system and solutions of the original system which remain near 
the stationary point. In the case of hyperbolic stationary points (defined below) 
this works very well. When a stationary point is not hyperbolic things are not 
so easy. In the latter case an important role is played by centre manifold theory. 
This theory will be illustrated by various cosmological examples in what follows. 
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One example of an application of centre manifold theory to homogeneous 
cosmological models in the literature can be found in |]]. There it was used to 
determine the asymptotic behaviour of spacetimes of Bianchi type VIIo with a 
radiation fluid. A feature of this example which is typical is that it represents a 
degenerate or borderline case. The case of a perfect fluid with equation of state 
p = (7 — l)p and 7 yf 4/3 can be analysed without centre manifold theory |l4j| 
and it is found that the behaviour changes qualitatively when passing through 
the case of the radiation fluid (7 = 4/3). 

In this paper two different applications of centre manifold theory are con¬ 
sidered. The first is to homogeneous spacetimes with collisionless matter. More 
specifically, it concerns locally rotationally symmetric (LRS) models of Bianchi 
type III. These play the role of a degenerate case within more general classes 
of Bianchi models with collisionless matter studied in |d| and 110 . The centre 
manifold analyses of the present paper extend and complete the results of j| 
and 0 - The other application is to homogeneous and isotropic spacetimes 
with a massive scalar field. The centre manifold analysis provides a rigorous 
confirmation of heuristic conclusions on inflation in this class of spacetimes due 
to Belinskii et. al. [||. 

The results on models with collisionless matter obtained in the following 
address a question raised in 0 - Suppose we have an expanding cosmological 
model with collisionless matter described by the Vlasov equation. The particles 
composing the matter have a certain non-zero velocity dispersion. It might be 
expected that if these particles are massive the velocity dispersion will decay 
with time as a consequence of the cosmological expansion so that the spacetime 
will look more and more like a dust model in the phase of unlimited expansion. 
It was shown in jn| that this is indeed the case for LRS models of Bianchi types 
I and II. For Bianchi type III, however, despite the amount of information on 
the dynamics obtained in JrJ, this issue was not settled. The reason for the 
difficulty is that in Bianchi type III one of the scale factors has an anomalously 
slow rate of expansion. In fact the results of 0 did not even suffice to show 
that this scale factor is unbounded. It is shown below that dust-like asymptotics 
also occurs in the case of Bianchi type III and that the problematic scale factor 
does increase without limit. An analogous result for massless particles is also 
obtained, thus completing the results of |J. It is to be hoped that these results 
will also open the way to further progress in understanding inhomogeneous 
models with collisionless matter such as those discussed in |JT] . 

The second section contains some general information about centre manifold 
theory. Some equations for spacetimes of Bianchi type III with kinetic matter 
are recalled in Section [| Centre manifold theory is applied to Bianchi III 
spacetimes with collisionless matter in the case of massive particles in Section 
[| and in the case of massless particles in Section |^. Section contains results 
on homogeneous and isotropic solutions with a massive scalar field. The final 
section presents some insights which can be obtained from the results in the 
body of the paper. 
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2 Centre manifold theory 

Let x = f(x) be a system of ordinary differential equations, where a; is a point 
of R” and the function / is smooth (C 00 ). Consider a point Xq with f(x o) = 0. 
Then x(t ) = xq is a time-independent solution of the system. In other words, 
xo is a stationary point of the given dynamical system. The linearization of 
the system about xq is the system y = Ay where the matrix A has entries 
Aj = dfi /dx l . If each eigenvalue A of A has the property that the real part 
of A is non-zero then the stationary point is called hyperbolic. In that case 
the theorem of Hartman and Grobman (Q, p. 120) shows that the linearized 
system is topologically equivalent to the original system in a neighbourhood 
of Xq. In other words, there is a homeomorphism of a neighbourhood of xq 
onto a neighbourhood of the origin in R” such that solutions of the non-linear 
system are mapped onto solutions of the linearized system. Note that this 
homeomorphism cannot in general be improved to a C 2 diffeomorphism. (See 
®» P- 127.) 

When the matrix A has eigenvalues which are purely imaginary and, in par¬ 
ticular, when it has zero as an eigenvalue things are more complicated than for 
hyperbolic stationary points. Recall that the generalized eigenspace associated 
to an eigenvalue A of A consists of all complex vectors x such that ( A—XI) k x = 0 
for some positive integer k. The centre subspace E c of the dynamical system 
at xq is defined as follows. Take the direct sum of the generalized eigenspaces 
corresponding to all purely imaginary eigenvalues and then restrict to the real 
subspace. Then E c is an invariant subspace of A and the restriction of A to E c 
has purely imaginary eigenvalues. E c is the maximal subspace having this prop¬ 
erty. The centre manifold theorem || says that there is a manifold M c passing 
through Xq which is invariant under the flow of differential equation and whose 
tangent space at Xq is E c . The manifold M c is called a centre manifold. It can 
be chosen to be C k for any positive integer k but not necessarily to be C°°. 
(See the example on p. 29 of ||.) It is also in general not unique. However any 
two centre manifolds of class C k are tangent to order k at xo and in fact the 
coefficients in the Taylor expansion of this submanifold at £o can be calculated 
recursively. Examples of this procedure will be seen later on. 

There is a generalization of the Hartman-Grobman theorem to non-hyperbolic 
stationary points called the reduction theorem J^| which reduces the study of 
the qualitative behaviour of solutions of a system near a stationary point to the 
study of solutions on the centre manifold. In favourable cases the qualitative 
behaviour of solutions on the centre manifold near the origin can be determined 
once a finite number of expansion coefficients of the centre manifold are known. 
This will be illustrated by the examples of cosmological models treated later. 

The reduction theorem by itself can provide useful information about the 
qualitative behaviour of solutions of a system of ordinary differential equations. 
For example this was done for the Mixmaster solution in ||. (Note, however, 
that since then a much deeper analysis of this system has been carried out 
in [|Ilj and |l2j.) This kind of straightforward application of the reduction 
theorem is not what is meant by ‘centre manifold theory’ in the present paper. 
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The latter should rather be reserved for the procedure where the existence of 
a centre manifold is used as a basis for doing calculations which give concrete 
information on the nature of solutions. It is calculations of this type which are 
the main tool of this paper. 

Centre manifold theory has the advantage of reducing the dimension of the 
dynamical system whose qualitative behaviour has to be determined. The di¬ 
mension of the centre manifold is of course just equal to the number of purely 
imaginary eigenvalues of the linearization, counting multiplicity. When the cen¬ 
tre manifold is one-dimensional things become particularly easy. Determining 
the stability of the origin in the centre manifold reduces to algebra. This hap¬ 
pens in the example studied in Section When the centre manifold is two- 
dimensional things are a lot more complicated but there is still a far-reaching 
theory available (see [7|, section 2.11) and in particular cases such as those stud¬ 
ied in Sections |] and 5 of this paper, direct approaches may lead to the desired 
result. 


3 The Bianchi III equations with massive parti¬ 
cles 

The results of this section extend the theorems of jl(| and the fundamental 
equations will be taken directly from that source. Here it will be assumed from 
the start that the Bianchi type is III. The spacetime metric is of the form 

ds 2 = -dt 2 + a 2 {t){9 1 ) 2 + b 2 {t){{9 2 ) 2 + ( 0 3 ) 2 ) , (1) 

where the 9 l are one-forms whose exterior derivatives satisfy the relations re¬ 
quired to ensure that the metric defines a spacetime of Bianchi type III. For 
instance we can assume the relations 99 1 = d9 2 = 0 and d9 3 = 9 2 A 9 3 . This 
metric is locally rotationally symmetric (LRS). The matter content of space- 
time is described by the spatially homogeneous phase space density of particles 
f(t,Vi). It is assumed that / is invariant under rotations in the (u 2 ,U 3 )-plane 
and the reflection v\ i—> —v\. The energy-momentum tensor Tij for the Einstein- 
Vlasov system with particles of mass m > 0 is diagonal and is described by 

p = J fo(.Vi)(m 2 +a~ 2 (v 1 ) 2 + b~ 2 ((v 2 ) 2 + {v 3 ) 2 )) 1/2 (ab 2 )~ 1 dvidv 2 dv 3 , 

Pi = J fo(vi)g n (vi) 2 (m 2 + a~ 2 (vi) 2 + b~ 2 {{v 2 ) 2 + (v 3 ) 2 ))~ l/2 (ab 2 )~ 1 dv 1 dv 2 d'ii^) 

where p is the energy density and pi = T l i the pressure components of the 
energy-momentum tensor. The function /o is determined at some fixed time 
to by fo(vi) = f(to,Vi) where / is the phase space density of particles. Some 
further technical conditions will be imposed on f 3 . It is assumed to be non¬ 
negative and have compact support. It is also assumed that the support does 
not intersect the coordinate planes Vi = 0. A function f 3 with this property 
will be said to have split support. This assumption ensures that the dynamical 
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system which describes the evolution of the spacetimes of interest has smooth 
coefficients. 

The momentum constraints are automatically satisfied for these models. 
Only the Hamiltonian constraint and the evolution equations are left. Instead 
of considering a set of second order equations in terms of a and b , we will refor¬ 
mulate these equations as a first order system of ODEs by introducing a new 
set of variables. The mean curvature tr k (where k l3 is the second fundamental 
form) is given by 

tr k = —( a~ 1 dci/dt + 2b~ 1 db/dt) . (3) 

A new dimensionless time coordinate r is defined by — 1 J* tr k(t)dt for some 
arbitrary fixed time <o- In the following a dot over a quantity denotes its deriva¬ 
tive with respect to r. The Hubble variable H is given by H = — trfc/3. This 
section and the next are concerned with the case of massive particles (m > 0). 
Define the following dimensionless variables: 

2 = m 2 /(a -2 + 26 -2 + m 2 ) , 

s = b 2 /(b 2 + 2 a 2 ) , 

M 3 = 36 -2 (trfc) -2 , 

S+ = —3(6 _ 1 d6/dt)(trfc ) _1 - 1 , (4) 

These variables lead to a decoupling of the equation for the only remaining 
dimensional variable H (or equivalently trfc) 

H = -(l + q)H, (5) 

where the deceleration parameter q is given by 

q = 2E^ + ±H(1 + R) . (6) 

The quantity R is defined by 

R= (pi + 2p 2 )/p , (7) 


where 

Pi Ip = (1 - z)sgi/h , 

Vi/p = |(1 - z)(l - s)g 2 /h , 

0i,2 = J fo{vi)(vi,2) 2 [z + {l-z)(s(v !) 2 + ±(l-s)((v 2 ) 2 + (v 3 ) 2 ))]- 1 / 2 dv 1 dv2dv 3 , 
h = J fo{vi)[z + (1 - 2 )(s(ul ) 2 + 1(1 - s)((u 2 ) 2 + (v 3 ) 2 ))] 1/2 dv 1 dv 2 dv 3 . ( 8 ) 

The assumption of split support ensures that the function R(s, z) is a smooth 
(■ C°°) function of its arguments. The related quantity R + defined by 

R+ = {P 2 ~Pi)/p ( 9 ) 
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is a smooth function of s and z for the same reason. 

The normalized energy density fl = p/(3H 2 ) is determined by the Hamilto¬ 
nian constraint and, in units where G = 1 /87r, is given by 

n = 1 - - m 3 . (10) 

The assumption of a distribution of massive particles with non-negative mass 
leads to inequalities for R, R + and fi. Firstly, 0 < R < 1 with R = 0 only when 
z = 1 and R = 1 only when z = 0. Secondly, —R < R + < ^R with R + = i R 
for s = 0 and R + = —R for s = 1. Thirdly f2 > 0. Using these inequalities in 
equation (j(|) in turn results in 0 < q < 2 (i-e., the same inequality as for causal 
perfect fluids, see Ju|). 

The remaining dimensionless coupled system is: 

S-j_ = —(2 — <?)^+ + -If3 T Ui?_1- , 
s = 6s(l — s)£+ , 
i = 2z(l — z)( 1 + E + — 3S + s) , 

M 3 = 2(«-S + )M 3> (11) 

In the following this is referred to for short as the Bianchi III system. It is of 
interest to note that the metric functions a, b are expressible in terms of s, z in 
the massive case. The relations are 

a 2 = z(m 2 s(l — z))~ x , b 2 = 2z(m 2 (l — s)(l — z))~ l . (12) 

There are a number of boundary submanifolds: 

z = 0,1 , 

s = 0,1 , 

0 = 0. (13) 

The submanifold z = 0 corresponds to the massless case. The submanifold 
z = 1 leads to a decoupling of the s-equation, leaving a system identical to 
the corresponding dust equations. The submanifolds s — 0, s = 1 correspond 
to problems with f 0 being a distribution while U = 0 constitutes the vacuum 
submanifold with test matter. 

Including these boundaries yields a compact state space. In order to apply 
the standard theory of dynamical systems the coefficients must be C 1 on the 
entire compact state space G of a given model. This is necessary even for 
uniqueness. In the present case it suffices to show that R and R+ are C 1 on 
G, i.e, that they are C 1 for s,z when 0<s<l,0<z<l. As has already 
been pointed out, this follows from the assumption of split support, which even 
implies the analogous statement with C 1 replaced by G°°. It would be possible 
to get C 1 regularity under the weaker assumption that /o vanishes as fast as 
a sufficiently high power of the distance to the coordinate planes but this is of 
little relevance to the main concerns of this paper. 
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Of key importance is the existence of a monotone function in the ‘massive’ 
interior part of the state space: 

M = (s(l-s) 2 )“ 1 / 3 2 (l- 2)- 1 , 

M = 2 M . (14) 

Note that the volume ab 2 is proportional to M 3 / 2 . This monotone function rules 
out any interior u>- and a-limit sets and forces these sets to lie on the s = 0 , 
s = l,z = 0 orz = l parts of the boundary. 

The physical state space, G , of the LRS type III models is given by the 
region in R 4 defined by the inequalities M 3 >0, 0<s<l,0<z<l, and 
1 — £+ — M 3 >0. A solution of the Bianchi III system (|Tl]) will be said to lie 
in the physical region if it lies in the interior of G. In the following we will be 
interested in the behaviour in the limit r —> oo of the solutions of the Bianchi III 
equations which lie in the physical region . This corresponds to the asymptotics 
of the underlying cosmological model in a phase of unlimited expansion. It 
was shown in that as r tends to infinity a solution in the physical region 
converges to a point of the line of stationary points L 3 which consists of all 
points of the form (^, 1, zq, |). The point to which the solution converges must 
satisfy zo > 0. It was left open in jllj whether that point must satisfy zo = 1 
or whether values less than one are also possible. In the next section it will be 
shown that in fact z —> 1 as r —> oo. 


4 Centre manifold analysis for the Bianchi III 
equations with massive particles 

It follows from Theorem 5.1 of []T(| that any solution of the Bianchi III system in 
the physical region converges to a limit of the form (^, 1 , zq, |) with 0 < z 0 < 1 
as r —> oo and that no solution in the physical region converges to a limit of 
this form as r —> — oo. In order to decide which values of z 0 can occur as 
limiting values of solutions in the physical region in the forward time direction 
it is necessary to study the behaviour of solutions close to the stationary points 
on the line L 3 in detail. If the system is linearized about one of these points 
then it is seen that the linearization has two negative eigenvalues and two zero 
eigenvalues. The fact that there is a whole line of stationary points implies that 
there is automatically one zero eigenvalue. However the second zero eigenvalue 
means that there is an extra degeneracy. For zq = 1 the matrix defining the 
linearized system is diagonalizable while for 0 < zq < 1 it is not. For this reason 
there are two significantly different cases to be considered. 

Consider first the case Zq = 1. New coordinates will be introduced in order 
to simplify the analysis. First the stationary point with zq = 1 will be translated 
to the origin. Define x = E + — y = 1 — s, w = M 3 — z = z — 1. Then the 
transformed system is 

x = — \ + q{\ + x) — 2x + w + HR + (15) 
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V= —6j/(l — y)(| + x) 

(16) 

z = —2(1 + z)z(—2x + | y + 3 xy) 

(17) 

w = 2(q-x~ i)(w + |) 

(18) 

q = ^ T 2x A 2x 2 A ^12(1 A R ) 

(19) 

Q = —x — w — x 2 

( 20 ) 


Next set Z = —z, u = x A w and v = x — w. In order to determine the 
behaviour of solutions of the resulting system near the origin we attempt to 
use the linearization at that point. The linearization has zero as an eigenvalue 
with a corresponding two-dimensional eigenspace. This eigenspace is the centre 
subspace. Centre manifold theory |dj tells us that there is an invariant manifold 
through the origin whose tangent space is the centre subspace. This manifold 
is a centre manifold. A centre manifold does not in general need to be unique. 
However approximations to it can be derived and these often suffice to determine 
the qualitative behaviour of solutions on the centre manifold. The qualitative 
behaviour of solutions near the stationary point is determined by the qualitative 
behaviour of solutions on any centre manifold. A general property of centre 
manifolds is that a stationary point has an open neighbourhood such that any 
other stationary point in that neighourhoocl lies on any centre manifold of the 
original point. It follows that in a neighbourhood of the stationary point of the 
system under consideration the part of the line L 3 close to that point lies on 
any centre manifold. In the case of present interest the centre manifold can be 
defined by the equations y = <p(u , Z) and v = ip(u, Z) where <f> and ip vanish at 
least quadratically at the origin. The invariance of the centre manifold implies 
that 


y = (d<p/du)u + (d<p/dZ)Z (21) 

v = ( dip/du)ii A ( dip/dZ)Z (22) 

on the centre manifold. The right hand sides of these equations vanish to at 
least third order at the origin. On the centre manifold 

y = —3 <p — 2><p{u + ip) + 3(p 2 (l + u + ip) (23) 

Substituting (|23|) into (||^) shows that (p vanishes at quadratic order. Putting 

this information back into the equation ( | 2 l| ) shows that </> vanishes at third 

order. This can be repeated indefinitely, with the result that (p vanishes to all 
orders. On the centre manifold 

= - §ip - | u 2 A \uip A | ip 2 - \u 3 A \u 2 ip + nip 2 + | ip 3 
+n(R + -±R) + ±n(l + R)(3ip-u) (24) 

Note that (R + R + )(y, Z) = 0(y) since R+R + is smooth and vanishes when y = 
0. Since (p vanishes to all orders it follows that R+ and —R are equal to all orders 
on the centre manifold. Hence to quadratic order ip(u, Z) = —j^u 2 — r 2 uZ +... 
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where r 2 is the derivative of R with respect to z, evaluated at the stationary 
point of interest. Next the dynamics on the centre manifold will be examined. 


u = u 2 + uip + \ip 2 + + i§«V + T6 u ^ 2 ~ fe^P 3 


+fl[(R + R + ) + \R{3u-iP)] (25) 

while 

Z = 2 uZ - 2 uZ 2 + 2Zip - 2 Z 2 ip - 3Z(1 - Z)<p{ 1 + u + ip) (26) 

Discarding terms of order greater than two gives the truncated system 

u = u 2 (27) 

Z = 2 uZ (28) 


For this system u~ 2 Z is a conserved quantity and so the qualitative behaviour 
of the solutions is easily determined. Is the qualitative behaviour of solutions 
of the full system on the centre manifold similar? The terms on the right hand 
side of the restriction of the evolution equations to the centre manifold written 
out explicitly above contain a factor of u. This is not an accident and in fact 
the full equations are of the form 


u = uf(u, Z) (29) 

Z = ug(u, Z) (30) 

for some functions / and g of any arbitrary finite degree of differentiability. This 
is because of the fact, already mentioned above, that the line L 3 of stationary 
points is locally contained in the centre manifold. Thus the restriction of the 
system to the centre manifold has a line of stationary points with u = 0. This 
means that the vector field defining the dynamical system vanishes for u = 0. 
The existence of the functions / and g then follows from Taylor’s theorem. The 
fact that L 3 lies on the centre manifold also implies that ip{u , Z) = uh(u , Z) for 
a function h of arbitrary finite differentiability. 

It is possible to get more precise information about the function / as follows. 
Recall that R + i?+ vanishes modulo the function <p. Since <p vanishes to all 
orders it follows that expressions containing a factor R + R + can be ignored in 
all perturbative calculations. We have 

D = —u — ju 2 — |uip — j'tp 2 (31) 

It follows that every term in the evolution equation for u on the centre manifold 
contains a factor of u 2 , either directly or via the fact that ip = uh. Hence 

f( 0 ,z) = 0 . 

For the purpose of qualitative analysis we can replace the above system by 
the system 


(32) 

(33) 
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u = f(u, Z) + ... 
Z' = g{u, Z) + ... 



which has the same integral curves away from the Z- axis. Only solutions with 
Z > 0 are physical. The Z -axis is an invariant manifold of the rescaled dynam¬ 
ical system. This new system has a hyperbolic stationary point at the origin 
which is a source. This makes it easy to determine the phase portrait for the 
system near the endpoint of L3 at z = 1. What we see is that there is an open 
set of initial data in the physical region such that the corresponding solutions 
approach the stationary point with z = 1 as r —> 00 and there is a neighbour¬ 
hood of the endpoint of L3 at 2 = 1 such that in that neighbourhood no solution 
in the physical region approaches a point of L3 with z < 1. 

Consider now the case 0 < Zo < 1. The analysis is similar to that in the 
case already considered, although the calculations are somewhat heavier. For 


Zq < 1 the linearization has a non-diagonal Jordan form. Define x = E + — 
y = 1 — s, w = M3 — z = z — Zo- Then the transformed system is 

x = — \ + q(\ + x) — 2x + w + flR + (34) 

y =-&y(l - y){\ + x) (35) 

5 = 2 (zo + z)( 1 — zo — z)(—2x + \y + 3a :y) (36) 

w = 2(q-x-±)(w + l) (37) 

? = | + 2a; + 2a; 2 -I- ^12(1 + R) (38) 

O = —x — w — x 2 (39) 


The linearization has zero as an eigenvalue with a corresponding two-dimensional 
generalized eigenspace. 

For calculational purposes it is convenient to choose coordinates which re¬ 
duce the linearization to a simple form. For this reason the following variables 
will be defined. They almost reduce the linearized operator to Jordan form. 


U = -2z 0 (l - z 0 )(x + w) (40) 

V = -l[(R-l)x + (R+l)w] (41) 

Y = y (42) 

Z = 5 + z 0 (l - z 0 )y +§[(-&- l)x + (R + l)tu]z 0 (l - z 0 ) (43) 

Here R is the value of R when s = 1 and z = Zq. Note that the value of i?+ 
when s = 1 and z = Zq is —R. Up to linear order the equations in the new 
coordinates are 

U = 0 + ... (44) 

u = -fu + ... (45) 

Y = -3 Y + ... (46) 

Z ={R+l)U + ... (47) 

In the new coordinates the centre subspace is defined by V = 0 and Y = 0. 
The centre manifold is defined by equations of the form Y = Z) and 


V = Z), where </> and ip are functions of any desired finite degree of dif¬ 

ferentiability which vanish at least quadratically at the origin. (These are not 
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the same functions <f> and ip introduced in the centre manifold calculation for 
zq = 1.) Now apply the equations 

Y = {d(f>/dU)U + (, d(j)/dZ)Z (48) 

V = (dip/dU)U + (dip/dZ)Z (49) 

which must hold on the centre manifold. The linear contribution to Z is pro¬ 
portional to U and hence the quadratic contribution to the right hand side of 
the first equation contains only terms proportional to U 2 and UZ and no term 
proportional to Z 2 . The quadratic part of the left hand side is proportional to 
the quadratic part of qb. Suppose that the latter is given by AU 2 + BUZ + CZ 2 
for constants A, B and C. Comparing the coefficients of Z 2 shows that C = 0. 
Then comparing coefficients of UZ shows that B = 0. At this point it has been 
shown that the right hand side of the equation vanishes to quadratic order and 
it follows that A = 0. Hence <j> vanishes to quadratic order. As in the case 
zo = 1 this information can be used repeatedly to show that <f> vanishes to all 
orders at the origin. 

Next the equation for V on the centre manifold will be examined. It turns 
out that the quadratic contribution to V is a linear combination of U 2 and UZ. 
It can be concluded that up to quadratic order ip = aU 2 + f3UZ + ... for some 
constants a and (3. These constants depend on the value of Zq at the point 
where the linearization is carried out and on the initial data /o- 

Substituting the information about the centre manifold into the original 
system shows that the restriction of the dynamical system to the centre manifold 
has the following form up to quadratic order: 

U = 7C/ 2 + ... (50) 

Z = {R+l)U+ SU 2 + eUZ + ... (51) 

for some constants 7, 6 and e. 

All the terms on the right hand side of the restriction of the evolution equa¬ 
tions to the centre manifold written out explicitly above contain a factor of U 
and for the same reason as in the case zo = 1 the full equations are of the form 

U = Uf(U,Z) (52) 

Z = Ug(U,Z ) (53) 

for some functions f and g of arbitrary finite differentiability. Thus the restric¬ 
tion of the system to the centre manifold has a line of stationary points with 
U = 0. 

Consider now the dynamical system defined by 

U' = f(U,Z) (54) 

Z' = g(U,Z) (55) 

where the prime denotes cL/dcr for a time coordinate cr. This has the same 
integral curves as the system we want to analyse in the complement of the Z- 
axis. It will be referred to in the following as the rescaled system. At the origin 
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/ vanishes while g has the non-vanishing value i?+ 1. In fact /(0, Z) = 0 for all 
Z, the argument being just as in the case zq = 1. As a consequence the integral 
curve of the rescaled system passing through the origin lies on the Z- axis. Since 
the vector field defining this system is non-vanishing at the origin it is clear that 
no solution can approach the origin from the physical region. 

Putting together the results of the two centre manifold analyses with Theo¬ 
rem 5.1 of yields the following theorem: 

Theorem 1 If a smooth non-vacuum reflection-symmetric LRS solution of 
Bianchi type III of the Einstein-Vlasov equations for massive particles is rep¬ 
resented as a solution of ( 0 ) then for r —> oo it converges to the point with 
coordinates (|, 1,1, |). The quantities pi/p and P 2 I p converge to zero as t —> 00. 

Next some information will be obtained on the detailed asymptotic behaviour 
of the spacetime as t —> 00. Consider first the behaviour of a solution of the 
rescaled dynamical system for zo = 1 as a —> — 00. At the origin the rescaled 
system has a stationary point which is a hyperbolic source. By a theorem of 
Hartman (see [;7j, p. 127) the system can be linearized by a C 1 mapping. As 
a consequence u = uoe a + o(e a ) for some constant uq and Z = o(e a ). Along a 
solution of the equations dr/da = u~ 1 (a). Hence for large negative values of a 
we have r = — u// 1 e~‘ T +o(e~ a ). Hence on the centre manifold u = — t _ 1 +o(t~ 1 ) 
and Z = o(t -1 ). 

General results on centre manifolds || show that the general solution satisfies 

u(t) =—t~ x + o^ -1 ) (56) 

z(t) = 1 + o(r _1 ) (57) 

y(r) = o(r~ k ) (58) 

V ( T ) =-J2 r ~ 2 +°( T ~ 2 ) ( 59 ) 

where k is any positive integer. As r —> 00 we have q = \ + o(l). Now 
dH/dt = — (1 + q)H 2 . It follows that H = |f -1 (l + o(l)). Thus in leading order 
r = | logt + .... Putting this information in the asymptotic expressions above 
gives the leading order behaviour of u, z, y and v as functions of t. This in turn 
gives the asymptotic behaviour of the variables E+, s, M 3 and 2 . The result is 


E+(t) = i^|(logt )- 1 + ... (60) 

s = 1 + o((logf) -fc ) (61) 

M 3 = | — |(logt ) _1 + ... (62) 

2 = 1 + o((logt) _1 ) (63) 


From this it can be concluded using (|li|) that a/ (logf ) 1 / 2 goes to infinity as 
t —> 00 . It follows directly from the fundamental equations that a goes to 
infinity slower than any power of t. To see this note that since tends to 
the quantity ( 6 - 1 d 6 /df)(trfc ) -1 tends to Hence b~ 1 db/dt = t~ l + ... 
and a~ 1 da/dt = o(t _1 ). This gives the desired conclusion concerning the upper 
bound for the rate of growth of a. The leading order behaviour of b can be read 
off from the defining equation for M 3 , with the result that b = t + .... 
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5 Centre manifold analysis for the Bianchi III 
equations with massless particles 

The analysis of LRS Bianchi spacetimes with massless particles in Pi included 
the case of Bianchi type III. In this section the results obtained in |J will be 
sharpened. As previously indicated, the dynamical system for LRS Bianchi type 
III spacetimes with massless particles is equivalent to the boundary component 
z = 0 oi the region G on which the system for massive particles is defined. 
Here we adopt the notation of [|lo| rather than that of [j|] . Then the dynamical 
system required in this section can be obtained by specializing a system we had 


in the last section. The result is 

x = —j + q(\ + x) — 2 x + w + f IR + (64) 

V = - 6 y(l - y)(\ + x) (65) 

w = 2 {q — x - 7 })(w + |) ( 66 ) 

q = \ + 2x + 2x 2 + fl (67) 

fl = —x — w — x 2 ( 68 ) 


The only modifications are that the equation for z has been dropped and that the 
relation R = 1 has been used. This is a three-dimensional dynamical system. Its 
linearization at the origin has eigenvalues 0, — | and —3 and its kernel is spanned 
by (1,0,0). It follows that this point has a one-dimensional centre manifold. 
It was already shown in || that all solutions of this system corresponding to 
solutions of the Einstein-Vlasov equations with massless particles converge to 
the origin as r —» oo. The purpose of the centre manifold analysis here is to 
obtain more detailed information about how this point is approached and hence 
how the scale factors behave in the phase of unlimited expansion. It is useful to 
substitute in the expressions for Q and q, thus obtaining the following dynamical 
system: 

x = \w + |* 2 — xw + x 3 — (x + w + x 2 )(l + R + ) (69) 

y = ~6y(l - y)(\ + x) (70) 

w = 2(-w + x 2 )(w + |) (71) 

The centre manifold is defined by the equations y = <p(x) and w = tj)(x). Us¬ 
ing the same procedure as in the last section the functions </> and ip can be 
determined up to quadratic order. In this case it turns out that <p vanishes to 
quadratic order while to that order ip = x 2 + .... Substituting this information 
into the evolution equation for x shows that x = 4x 2 + 0(|x| 3 ) on the centre 
manifold. Since physical solutions converge to the origin as r —> oo they must 
be approximated by solutions on the centre manifold with x negative. It can be 
concluded that for these solutions E + is eventually less than one half. 

With the information already obtained, various features of the dynamics 
of the spacetime can be reconstructed. To leading order x = —(1/4r) + ..., 
y = o(t~ 2 ) and w = l/16r 2 + o(r -2 ). The quantity q tends to \ and so some 
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of the calculations for the massive case can be taken over here. As in that case 
b = t + .... It is not possible to determine the asymptotic behaviour of the 
scale factor a in the same way as in the massive case; a different approach is 
necessary. It may be computed that a~ 1 da/dt = ^(trk)x. Hence in leading 
order a~ 1 da/dt = ^(tlogf ) -1 + .... It follows that a —> oo as t —> oo. 

6 Inflation with a massive scalar field 

This section presents a rigorous derivation of some heuristic results of [|| for 
spatially flat homogeneous and isotropic spacetimes with a massive scalar field 
as source. The starting point is the following system of equations from [|j. 


x v = y (72) 

y v = ~x- 3 y(x 2 + y 2 ) 1/2 (73) 

This is a dynamical system on the plane whose coefficients are smooth every¬ 
where except at the origin, where they are C 1 . Transforming to polar coordi¬ 
nates (r, 9) and introducing p = r/(l +r) gives the system 

p v = —3 p 2 sin 2 9 (74) 

9 V = — 1 — 3p(l — p ) -1 sin 0 cos 0 (75) 

Introducing a new time coordinate r and setting u = 1 — p gives 

u T = 3m(1 — u) 2 sin 2 9 (76) 

0 T = —u — 3(1 — u) sin 0 cos 0 (77) 


So far this is nothing new compared to what is done in 0. Now the stationary 
point at the origin will be investigated. The eigenvalues of the linearization 
are —3 and zero. Evidently the 0-axis is invariant and in fact it is the stable 
manifold of the origin. The centre subspace is given by it + 30 = 0 and thus we 
introduce v = u + 30 in order to study the centre manifold. This leads to the 
transformed system 


u T = 3it(l — u ) 2 sin 2 (|(i; — it)) (78) 

v T = — 3mcos 2 (|(w — u)) — 6 u 2 sin 2 (i(t; — it)) + 3m 3 sin 2 (|(i; — u )) 

—9sin(|(ti — it)) cos(|(t; — it)) -f 9wsin(l(i; — it)) cos(|(u — it)) (79) 

The centre manifold is of the form v = 4>[u). To quadratic order = —it 2 + 
.... On the centre manifold u T = ^it 3 + 0(it 4 ) and hence there is a unique 
solution which enters the physical region. 

The asymptotic form of the solution asr-» —oo is it = a/3/2(—t) -1 / 2 +.... 
It follows that v = + .... Hence 

P = 1 - VsJ^i-T )- 1 / 2 + ... (80) 
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and 


0 = -v/i 76 (-t )- 1/2 + ... (81) 

The next step is to convert back to the original variables. In leading order 

r] = — V&(— r) 1 / 2 + ... and so p = 1 + 3/rj + ... and 0 = 1 /677 +_ It follows 

that x = 77/3 + ... and y = jg + .... The combination z = \Jx 1 + y 2 occurring 
in J|] is asymptotic to 77/ 3 . The linear dependence of x and 2 with respect to 77 
is what was found in Q. The interpretation of these variables is that 77, x and 
0 are proportional to proper time, the scalar field </> and the mean curvature of 
the hypersurfaces of constant time respectively. Thus in the limit t —> — 00 both 
(f> and the mean curvature are proportional to t. The leading order behaviour 
of the scale factor for large negative times is, up to inessential constants, e* . 
This gives a rigorous confirmation of the conclusions of ||. Note that the 
relevance of centre manifolds in this context has been pointed out by Foster 
[[| but that he did not carry out a full centre manifold analysis; his paper was 
mainly concerned with other applications of dynamical systems to homogeneous 
and isotropic solutions of the Einstein equations coupled to a scalar field with 
potential. 


7 Further comments 

This section contains some further comments on the results of Sections [| and |j. 
In terms of the dimensionless variables which are the unknowns in the dynami¬ 
cal system studied in those sections, solutions of the equations with collisionless 
matter (with both massless and massive particles) converge at late times to a 
point corresponding to a vacuum solution of Bianchi type III. The corresponding 
statement holds for LRS Bianchi type III solutions with dust. (Cf. the discus¬ 
sion in section 6 of |J.) In particular, the dimensionless density parameter fl 
tends to zero in the expanding direction. Furthermore, in the case of massive 
particles, the ratios T ?/p of the spatial eigenvalues of the energy-momentum 
tensor (i.e. the principal pressures) to the energy density tend to zero. Thus 
in a certain sense both solutions with collisionless matter and dust solutions 
are approximated by vacuum solutions, while solutions with collisionless matter 
and massive particles are approximated better by dust solutions than either of 
these are approximated by vacuum solutions. This is one of the main results of 
this paper. 

One of the scale factors, a, grows much more slowly than the other scale 
factor b. An important conclusion is that a tends to infinity at all in the case 
of collisionless matter, both in the massive and massless cases. Centre mani¬ 
fold theory was the essential tool which allowed us to prove this and thus to 
go beyond what was achieved in || and jTo| . In this sense both types of colli¬ 
sionless matter resemble each other and differ from the vacuum case, where a 
is asymptotically constant. These results lend further support to the following 
two suggestions made in Q and The first is that any spatially homoge¬ 
neous solution of the Einstein equations with collisionless matter and massive 
particles behaves asymptotically in a phase of unlimited expansion like a dust 
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model. The second is that in any non-vacuum spatially homogeneous model 
with perfect fluid or collisionless matter eventually every scale factor increases 
in the time direction in which the volume increases. It would be interesting to 
investigate the truth of these statements in solutions of the Einstein equations 
with collisionless matter which are LRS Bianchi type VIII or type II but not 
LRS. 
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